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Takashi Hotta^ and Hisatomo Harima^ 

^Advanced Science Research Center, Japan Atomic Energy Agency, Tokai, Ibaraki 319-1195 
^Department of Physics, Faculty of Science, Kobe University, Kobe 657-8501 

(Received February 6, 2008) 

We propose an effective model on the basis of a j-j coupling scheme to describe local /-electron 
states for realistic values of Coulomb interaction U and spin-orbit coupling A, for future development of 
microscopic theory of magnetism and superconductivity in /"-electron systems, where n is the number of 
local f electrons. The effective model is systematically constructed by including the effect of a crystalline 
electric field (CEF) potential in the perturbation expansion in terms of 1/A. In this paper, we collect all the 
terms up to the first order of 1/A. Solving the effective model, we show the results of the CEF states for 
each case of n=2^5 with Oh symmetry in comparison with those of the Stevens Hamiltonian for the weak 
CEF. In particular, we carefully discuss the CEF energy levels in an intermediate coupling region with 
A/[/ in the order of 0.1 corresponding to actual /-electron materials between the LS and j-j coupling 
schemes. Note that the relevant energy scale of U is the Hund's rule interaction. It is found that the CEF 
energy levels in the intermediate coupling region can be quantitatively reproduced by our modified j-j 
coupling scheme, when we correctly take into account the corrections in the order of 1/ A in addition to the 
CEF terms and Coulomb interactions which remain in the limit of A=oo. As an application of the modified 
j-j coupling scheme, we discuss the CEF energy levels of filled skutterudites with Th symmetry. 

KEYWORDS: Crystalline electric field, spin-orbit interaction, j-j coupling scheme, LS" coupling scheme, 
filled skutterudites 



1. Introduction 

Theory of crystalline electric field (CEF) has been devel- 
oped so far on the basis of an LS coupling scheme. '"^ Nowa- 
days the CEF analysis is reduced to an automatic procedure, 
since we can simply refer the Hutchings table for CEF pa- 
rameters,^ B^, depending on the value of total angular mo- 
mentum J of relevant rare-earth or actinide ion. For instance, 
for Ce^+ ion with one / electron, due to a spin-orbit coupling, 
the ground state is included in the J=5/2 sextet. By consult- 
ing the Hutchings table for J=5/2,^ we easily find the matrix 
elements for CEF potential among six states of J =5 12. For a 
cubic system with Oh symmetry, after some algebraic calcu- 
lations, we immediately understand that the J=5/2 sextet is 
split into doublet and quartet. 

For the purpose to fit the experimental results on /-electron 
materials in a high-temperature region, the CEF analysis us- 
ing the Hutchings table is quite convenient, since it is not nec- 
essary to have deep knowledge on the origin of B'^, which is 
given by the sum of electrostatic potentials from the ligand 
anions surrounding the rare-earth or actinide ion. In fact, for 
the case of Pr'^+ ion including two / electrons, it is enough 
to consider the CEF effect for ,7=4 nontet, which is obtained 
after the consideration of the Hund's rules and the spin-orbit 
interaction in the LS coupling scheme. Thus, the competition 
between CEF effect and Coulomb interactions does not seem 
to appear in the actual calculation. This fact reduces our task 
in the CEF analysis. 

If / electrons are perfectly localized, we can satisfy with 
the above treatment and further improvement on the CEF the- 
ory may not be required. However, when /-electron proper- 
ties are gradually changed from localized to itinerant nature, 
there always appear rich phenomena including exotic mag- 
netism and unconventional superconductivity in /-electron 
systems. In order to describe the low-energy excitation of / 
electrons, we should improve the CEF theory so as to be com- 



patible with the /-electron Bloch state. For such a purpose, 
the LS coupling scheme is inconvenient, since the CEF po- 
tential is applied to the multi- /-electron state which is firmly 
constructed from the Hund's rules and the spin-orbit interac- 
tion. Rather, a j-j coupling scheme is more convenient, since 
individual /-electron states are first defined under the effect 
of CEF potential. Then, we construct the multi-/-electron 
state by considering the effect of Coulomb interactions, using 
standard quantum-field theoretical techniques. In contrast, in 
the LS coupling scheme, we cannot use such standard tech- 
niques, since Wick's theorem does not hold. 

In order to describe the /-electron system by the Bloch 
state, it is standard to exploit the relativistic electron band the- 
ory, which is based on the j-j coupling scheme. In this sense, 
it is quite natural to construct an /-electron model on the basis 
of the j-j coupling scheme. By keeping the same /-electron 
basis in a simple tight-binding model as that of the relativis- 
tic band-structure calculation, we can determine the hopping 
amplitudes in the tight-binding model to reproduce the Fermi- 
surface structure of the relativistic band-structure calculation 
result. After that, we attempt to include electron correlations 
between / electrons in the tight-binding model. 

On the above background, it has been recently proposed to 
construct a microscopic model for /-electron compounds on 
the basis of the j-j coupling scheme,"^' ^ by applying a tight- 
binding approximation for kinetic part of / electrons. Then, 
microscopic theories have been developed for the understand- 
ing of novel magnetism^'' and unconventional superconduc- 
tivity'""'^ in /-electron systems. In addition, it is also possible 
to study complex multipole phenomena from a microscopic 
viewpoint by using the same /-electron model. '^^' 

However, in the standard j-j coupling scheme in which 
j=7/2 states are simply discarded, only second- and fourth- 
order CEF parameters (S", B^, and are included. Be- 
cause of the symmetry reason, the effect of sixth-order terms 
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(-Bg and -Bg) are dropped. Thus, doublet does not ap- 
pear as a stable ground state in the j-j coupUng scheme for 
Oh symmetry. Moreover, we cannot include the effect of -Bg 
which is characteristic of filled skutterudites with Th symme- 
try.^^ These points should be improved for the further devel- 
opment of the microscopic theory of /-electron systems on 
the basis of the j-j coupling scheme. 

In this paper, we propose an effective model on the ba- 
sis of the j-j coupling scheme, which describes low-energy 
/-electron states by considering the effect of CEF potentials 
within the first order of 1/A, where A is the spin-orbit interac- 
tion. Since the effect of sixth-order CEF terms (J5§ and 
are included as two-body potentials for / electrons, we can 
reproduce the CEF energy levels in the intermediate coupUng 
region on the basis of the concept of the j-j coupling scheme. 

The organization of this paper is as follows. In Sec. 2, first 
we show the local /-electron Hamiltonian H composed of 
the spin-orbit couphng, the CEF potential, and the Coulomb 
interactions. On the basis of H, we discuss the /-electron 
states in the LS and j-j coupling schemes. Then, we derive 
the effective model from H by the perturbation expansion in 
terms of 1/A. In order to compare the results of the effec- 
tive model with those of the weak CEF region, we also define 
the Stevens Hamiltonian by the method of operator equiva- 
lents.^ In Sec. 3, we show the results of the effective model for 
the cases of n=2~5 in comparison with those of the Stevens 
Hamiltonian, where n denotes the local /-electron number. 
For each n, we compare the CEF energy levels of Hes with 
those of H in the intermediate coupling region. We also show 
that the CEF energy levels for Th symmetry are reproduced in 
our modified j-j coupling scheme. In Sec. 4, we briefly dis- 
cuss the /-electron state for n>7 in the j-j coupUng scheme, 
by showing the CEF energy levels for the cases of n=12 and 
13. We remark a couple of future issues concerning the appli- 
cation of our modified j-j coupling schemes. Finally, the pa- 
per is summarized. Throughout this paper, we use such units 
as ?i=A;b=1 and the energy unit is eV. 

2. Formulation 

2.1 Hamiltonian 

In general, the local /-electron Hamiltonian is given by 

H = Hso + -ffcEF + Hint- (1) 

The first term indicates the spin-orbit coupling, written as 

where A is the spin-orbit interaction, fma is tiie annihilation 
operator of / electron, (t=+1 (—1) for up (down) spin, m is 
the ^-component of angular momentum £=3, and the matrix 
elements are expressed by 

Cm,a;m,cr — ?TlCr/2, 

Cm+a,-a;m,a- = \/ ^(^ + 1) " + Cr)/2, 

and zero for other cases. 

The second term denotes the CEF potential, given by 

-f^CEF = ^ Bjn,m'fm<Tfm'a, (4) 

where Bm,m' is determined from the CEF table for J=£=3} 
Note that electrostatic CEF potentials do not act on /-electron 
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spin. For the cubic system with Oh symmetry, Bm.m' is ex- 
pressed by two CEF parameters for J=3, and Bg, as 

^3,3 = S_3,-3 = 180^2 + 180^0, 
B2,2 = B_2.-2 = - lOSOBg, 

Bi,i = = 605° + 27005°, 

Bofi = 360^2 - 3600BJ?, 
£3!-! = -B-3,1 = 60VT5(B0 - 21^0), 
52,-2 = 300^2 + 7560Bg. 

Note the relation of Bjn,m'=Bm' ,m- Following the traditional 
notation,^ we define 

El = Wx/F{A), Bl = W{1 - |x|)/F(6), (6) 

where x and the sign of W specify the CEF energy scheme, 
while \ W\ determines the energy scale for the CEF potential. 
Concerning non-dimensional parameters, F(A) and F(6), we 
choose F(4)=15 and F(6)=180 for ,7=3 by following Ref. 2 
In actual /-electron materials, the magnitude of the CEF po- 
tential is considered to be 10~^~10~^eV. For the calculation 
of the CEF energy level in this paper, we set 14^=— 10~^eV, 
even if we do not explicitly mention the value of W. 

Finally, ilint denotes Coulomb interaction term, given by 

Hint— ^ ^ ^^-^mim2,m3m4fmiafm2cr'fm3a'fm4Ti (7) 

where the Coulomb integral Imim2,m3m4, is expressed by 

6 

= '^F''ck{mi,m4)ck{m2,m3). (8) 
fe=o 

Here F'' is the Slater-Condon parameter^^'^'* and Ck is the 
Gaunt coefficient^^' which is tabulated in the standard text- 
books of quantum mechanics. Note that the sum is Hmited 
by the Wigner-Eckart theorem to k=0, 2, 4, and 6. 

In principle, the Slater-Condon parameters and the spin- 
orbit interaction are determined so as to reproduce the spec- 
tra of rare-earth and actinide ions for each value of n, but in 
any case, the Slater-Condon parameters are considered to be 
distributed between 1 eV and 10 eV, while A is in the order of 
0. 1 eV. In this paper, in order to discuss conveniently the com- 
petition between the Coulomb interactions and the spin-orbit 
coupling, we parameterize F'^ as 

F° = lOU, F^ = 5U, F'^ = W, F^ = U, (9) 

where U is introduced as a scaling parameter. As long as we 
set the Slater-Condon parameters between 1 eV and 10 eV, 
the results obtained in the paper do not change qualitatively, 
even if we use different values for F*^. As for the meaning of 
U, roughly speaking, it indicates the order of the Hund's rule 
interaction Jh among / orbitals, which is considered to be in 
the order of a few eV. However, for the purpose to explain the 
change between the LS and j-j coupling schemes, we treat 
U asa free parameter in the following arguments. 

2.2 LS versus j-j coupling schemes 

In order to describe the low-energy /"-electron states of 
H, there are two typical approaches, the LS coupling and j-j 
coupling schemes. In the LS coupling scheme, first the to- 
tal spin S and total angular momentum L are formed by fol- 
lowing Hund's rules. After forming S and L, we include the 
effect of spin-orbit interaction XlsL-S, where \LS=^/n for 
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Fig. 1. Overlap integrals (a) and (b) (vI'I^I'ls) foi'"=l~13 for 

several parameter sets of U and A. 



n<7 and \ls=~^/ (14 — n) for n>7. The ground state mul- 
tiplet is specified by the total angular momentum J, given by 
J=L+S. From simple algebra, the ground-state multiplet is 
characterized by J=\L—S\ for n<7, while J=L+S for n>7. 

On the other hand, in the j-j coupling scheme, first we in- 
clude the spin-orbit coupling so as to define the one-electron 
state labelled by the total angular momentum j, given by 
j=s+£, where s denotes one-electron spin. For / orbitals with 
£=3, we immediately obtain an octet with j=7/2(=3H-l/2) and 
a sextet with j=5/2(=3— 1/2), which are well separated by the 
spin-orbit interaction. Note here that the level for the octet 
is higher than that of the sextet. Then, we take into account 
the effect of Coulomb interactions to accommodate n elec- 
trons among the sextet and/or octet, leading to the ground- 
state level in the j-j coupling scheme. For n=1^6, we ac- 
commodate / electrons in the j=5/2 sextet to construct the 
multiplet specified by J, while for n=7'^13, / electrons are 
accommodated in the j=7/2 octet. 

As intuitively understood from the above brief explanation, 
the LS and j-j coupling schemes become exact in the limit of 
U=oo and A=oo, respectively. In actual /-electron materials, 
as mentioned above, J7(~ Jh) is as large as a few eV. On the 
other hand, for rare-earth atoms, A is about 0. 1^0.2 eV, while 
it is 0.2^0.3 eV for actinide atoms. Thus, the actual situation 
is characterized by X/U in the order of 0.1. This intermedi- 
ate coupling region is close to neither the LS coupling limit 
{U=oo) and the j-j coupling limit (A=oo), but due to the re- 
lation of U>X in actual materials, the LS coupling scheme 
was widely used for the description of the local /-electron 
state. Contrary to such a historical trend, we believe that it is 
meaningful to reexamine the significance of the j-j coupling 
scheme and to provide an alternative systematic framework to 
understand the CEF state of /"-electron compounds. 

For the purpose, let us see how the wavefunction of H is 
changed by U and A. Here we suppress the CEF potential to 
focus on the competition between the Coulomb interaction 
and the spin-orbit coupling. In Fig. 1(a), we depict the over- 
lap integral for n=1^13 for several sets of A and 
U, where \'^) is the ground state of Hso + i?int and 
denotes the state at the limit of A=oo. Here we use A=10^ and 
U=l to express virtually the j-j coupling limit. The solid cir- 
cles denote the overlap integral between the LS and j-j cou- 
pling schemes, {'i>Ls\'^j-j), where I^'ls) denotes the state 
at the LS coupling limit. Note that we set A=l and U=10^, 
where U is large enough to arrive at the LS coupling limit. 
We find that {'i/Ls\'^j-j) exhibits a "V"-shaped function of 
n. It is natural, since the Hund's rule interaction works most 



effectively at half filling (n=7), which is the worst case for the 
j-j coupling scheme. 

In the intermediate coupling region for A=0.1~0.3 and 
U=l, the overlap integrals are increased with the increase of 
A. In particular, for n<7 and A=0.1, the values of {'i>\'i>j^j) 
are larger than 0.6. This lower limit becomes 0.8 and 0.9 for 
A=0.2 and 0.3, respectively. These facts suggest that the j-j 
coupling scheme is the good starting approximation to con- 
sider the /"-electron state for n<7. 

Note that the values of for n=5 and 6 become 

larger than that for n=4, leading to a shallow minimum around 
at n=3~4, when A is increased. It is intuitively understood by 
the electron-hole relation among the j=5/2 sextet in the j-j 
coupling scheme. For n=6, the singlet (J=0) of fully occu- 
pied j=5/2 state is considered to be very stable, since it is the 
closed-shell structure in the j-j coupUng scheme. For n=5, 
when we accommodate five / electrons in the j=5/2 sextet, it 
corresponds to the one-hole case. In the j-j coupling limit, we 
expect that the wavefunction for one-hole state is the same as 
that of one-electron case due to the electron-hole conversion. 
Thus, the overlap integral is rather increased for n=5 and 6, in 
comparison with the case of n=4, when we increase A. This 
fact also supports that we can construct the ground state in 
the intermediate coupling region on the basis of the j-j cou- 
pling scheme. Note that for n>7, the overlap integral is also 
increased with the increase of A and n. In the region of n>8, 
the values seem to approach unity gradually, but for n=7, it is 
increased only slowly. For the half-filling case, it is difficult 
to recommend the usage of the j-j coupling scheme. 

In Fig. 1(b), we show the overlap integral with the LS 
coupling wavefunction, defined by (^'I^'ls)^ for n=l~13 for 
several sets of A and U. The solid circles are the same as those 
in Fig. 1(a), since they denote {'i'j^jl'^Ls)- For A=0.1 and 
U=l, the ground state of H is well approximated by the LS 
coupling state except for n=6, since we find (^'|^'j_j)>0.9 
for n=l~5 and 7^13. When we increase A, the overlap inte- 
grals are totally suppressed. For A=0.3, it seems to be difficult 
to use the LS coupling scheme, in particular, for the region 
near n=6. 

In order to construct an effective model in the intermediate 
coupling region with X/U in the order of 0.1, it is necessary 
to improve one of the LS and j-j coupling schemes. We pre- 
fer one of the schemes, depending on the nature of the current 
problem. For n<6 and ?t.>8, we believe that the j-j coupling 
scheme provides us the good starting point to approach the in- 
termediate coupling region. As mentioned in the introduction, 
it is difficult to include the itinerant nature of / electrons in 
the LS coupling scheme. In order to describe the Bloch state 
of / electron, we should prefer the j-j coupling scheme for 
the description of the /-electron state. Then, we consider the 
local /-electron problem to construct the effective model by 
the perturbation expansion in terms of 1/ A from the limit of 

A=CX3. 

In the j-j coupling scheme, even if we include the effect of 
the order of 1/ A, the local /-electron state is still composed 
of the j=5/2 sextet. When the periodic system is discussed by 
using the effective model constructed from the j=5/2 states, 
we consider the hybridization between j=5/2 states and con- 
duction bands, although in actuality, there should also exist 
the hybridization of j=7/2 states. However, in the relativis- 
tic band-structure calculations for light rare-earth materials. 
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the j=5/2 sextet is occupied, while the j=7/2 octet is unoc- 
cupied, suggesting that / electrons in the j=5/2 sextet mainly 
contribute to the ground and the low-energy excited states. 
In fact, the Fermi-surface sheets are composed of j=5/2 elec- 
trons. Since we are interested in the electronic structure near 
the Fermi level, it seems natural to exploit the j=5/2 states for 
the construction of the effective model to discuss low-energy 
electronic properties of /-electron materials. 

2.3 Effective model due to the expansion in terms ofl/X 

Now let us construct the effective Hamiltonian in the j-j 
coupling scheme by including the effect of the CEP potentials 
in the perturbation expansion in terms of 1/A. Note that the 
energy scale |T4^| of the CEP potential is much smaller than 
U and A. As already mentioned, we set T4^=— 10^^ eV in this 
paper to be consistent with such a situation of the weak CEP. 

First we transform the /-electron basis between (m, a) and 
(j, jj) representations, connected by Clebsch-Gordan coeffi- 
cients, where j is the total angular momentum and ji is the 
^-component of j. Hereafter we use symbols "a" and "6" for 
j=5/2 and 7/2, respectively. When we define /j^ as the anni- 
hilation operator for / electron labeled by j and /x, the trans- 
formation is defined as 



(10) 



where the Clebsch-Gordan coefficient Cj^^-rn,a is give by 



Ca,ii;ti-a/2,cr = -tT\/(7/2 - g/x)/7, 



Cu 



V(7/2 + a/x)/7, 



^6,/i;/i— cr/2,(T 

and other components are zero. 

After the transformation, we express H as 



H = iJs( 



(11) 



(12) 



where the spin-orbit coupling term Hgo is diagonahzed as 

H.o = Y.^iflfj^, (13) 

with Aa=— 2A and A6=(3/2)A. The remaining part includes 
the CEF and Coulomb interaction terms as 

H' = -ffcEF + -ffinti (14) 

where .ffcEF and ifjnt are given by 
and 

■ffint = ^ ^ -^iil''/r2'f^3A»4/li/il/l2;i2/?3;'3/74/i4 5 (16) 

respectively. Note that B and I are the CEF potential and 
Coulomb interactions expressed in the basis of j and /i. 

In order to obtain the effective model, we use the degen- 
erate perturbation theory by treating H' as a perturbation to 
Hso. The effective model is written as 



(17) 



where H^^^ denotes the fcth-order term with respect to 1/A 

for the multiplet labeled by j. The zeroth-order term Ha 
indicates the model in the standard j-j coupling scheme."^' ^ 



(0) 



When we simply ignore all the terms including the symbol h 
(j=7/2), we obtain Hi°^ as 

Zj(0) _ rr(0) „(0) 
-"o — ^CEF ^ ^int • 

The CEF term is given by 

ff(0) _ Da,a j-f J- 

-'-'CEF ~ / J -"ai,a2JaaiJ««2 5 



(18) 



(19) 



where B^'_^^^ denotes the CEP potential for J=5/2, given by 

^±5/2,±5/2 = 60B4"(l,5/2), 

2). 

(20) 



^±3A±3/2 = -180i?4°(l,5/2), 
%V±i/2 = 12050(1,5/2), 
B"'",. = B^3/2,±5/2 = 60x/5B2(l,5/2). 



±5/2,=f3/2 

Here (n, J) denotes the CEF parameter for n and J, where 
n is the local /-electron number and J denotes the total angu- 
lar momentum of ground state multiplet. For n=l and J=5/2, 
B2il, 5/2) is related with B^ for J=£=3 in eq. (6) as 

S0(l,5/2) = (11/7)5°. (21) 

The meaning of the coefficient 11/7 will be explained later. 
The Coulomb interaction term is given by 

H^°^ = faa,aa j-t j-f f f (22) 

ai~a4 

where the Coulomb integral is expressed by three 

Racah parameters, E^, Ei, and E2,'^'^ which are related to 
the Slater-Condon parameters as 

FO - (80/1225)^2 - (12/441)^4, 



Eq 

El = (120/1225)^2 _^ (18/441)^4^ 
E2 = (12/1225)^^2 _ (i/44i)i?'4_ 



(23) 



Among Coulomb interactions, the Hund's rule coupling is the 
relevant energy scale for the construction of multi-/-state. In 
this case, among three Racah parameters, E2 plays a role of 
the Hund's rule coupUng in the j-j coupling scheme. Con- 
cerning the magnitude, by combining eqs. (9) and (23), we 
obtain £;2=(31/15)[//49=0.04[/ in the present parameteriza- 
tion. Note that the magnitude of the Hund's rule interaction 
in the j-j coupling scheme is generally reduced. In order 
to understand this point intuitively, it is convenient to con- 
sider a simple Hund's rule term, expressed as — Jhs^, where 
s denotes the spin of / electron. Here we note the relation 
s={gj — l)j, where gj is the Lande's g-factor and j indi- 
cates the total angular momentum. Since gj=6ll for j=5/2, 
we obtain s=— (1 /7)j. Namely, the original Hund's rule term 
is rewritten as — ( Jh/49)j'^ in the j-j coupling scheme. Thus, 
the Hund's rule interaction in the j-j coupUng scheme is re- 
duced as Jh/49 from the original value. As discussed later, 
this reduction has an important meaning to understand the 
rather wide applicabiUty of the effective model. 

Now we consider the correction term in the order of 1/ A, 
which is written as 



— -"CEF+-"inf 



The CEP term is given by 

r(i) _ V 

CEF — 2^ ^ai,o 



H, 



if aot\ f 0- 



(24) 



(25) 
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with 



■"ai ,a2 



^ {ai\ HcEF\P) {P\HcEF\a2) 

/3 



Aa — Af, 



(26) 



where |a)=/t^|0) and ia^z-j^^^ 
some algebraic calculations, we obtain 



=fL-j\0). For Oh symmetry, after 



-°.V2,5/2 



B 



(1) 



^(1) 

3/2,3/2 
^(1) 
-°l/2,l/2 

^5/2,-3/2 ^3/2,-5/2 



5/2,-5/2 = 5£8/6 + e7/6, 
-3/2,-3/2 = £8/6 + 5£7/6, 



(1) 



-1/2,-1/2 



£8, 



(27) 



= V5{ss - £7)76, 
where £7 and £3 are, respectively, given by 

e7 = -(^)'^(554° + 1265°f, 

and 



(28) 



(29) 



These values indicate the energy corrections for and Fg 

states at n=l, respectively. Note that such corrections are in 
the order of W^/X, as easily understood from the definition 

01 jD(xi^a2 • 

Concerning the 1/A-correction to the two-body potential, 

^hit expressed as 



ai~a4 



j(i) 

aia2iO;3a4 



/aai /aa2 f f*^3 / a 



(30) 



where the two-body potential is formally given by 



7(1) 

a!ia2, 0:304 



{aia2\H'\al3){(3a\H'\a3a4) 

Aa - Ab 



+ 



E 

V- (aia2|g'|/3i/32)(/32/?i|g'|a3Q4) (31) 



,/32 



2(A„ - Ab) 



- (aia2|-ffcEFlQ^3a4), 

where \a,a2)= /L./lajO), |a/3>= /i„/b^^|0), and |/3i/32)= 
/b/3i /fe/32 1'^^' th^t the final term in eq. (31) is needed to 
avoid the double-count of the contributions from Hfl^^, when 
we diagonalize Hcs for n>2. As we will see later, the two- 
body effective potential /^^^ plays a crucial role in the inter- 
mediate coupling region. It may be possible to obtain the an- 
alytic form of l'-^^ after lengthy and tedious algebraic calcu- 
lations, but in this paper, we evaluate l'^' only numerically, 
since our purpose here is to show that Hes actually works. 
The derivation of more convenient analytic form of I^^^ is 
one of future tasks. 

Now we discuss the range of the value of A, in which Hes 
works. First we note that i?cEF '■^^ order of W'^/X. Since 
this term is small in the order of |VF|/A in comparison with 
iJ^gp, it does not play important roles. The first two terms in 
eq. (31) also include the contributions in the order of W'^/X, 
but we find that they are exactly cancelled by the last term. 
Thus, the lowest contribution to the CEF potential from H^^^ 
is in the order of |W^|?7/A. Note that the terms in the order of 
?7^/A are not important, since they contribute to the energy 
shift of the ground-state multiplet. 



From the mathematical viewpoint of the convergence of 
the perturbation expansion, |T4^|C//A is thought to be smaller 
than \ W\, which is the energy scale of -f^cEF- Namely, we ob- 
tain the condition of U/X<^1, but under this condition, HeS 
cannot be used for realistic systems at the first glance, since 
U is larger than A in /-electron compounds. In order to re- 
consider this point, we note that the sixth-order contributions 
of the CEF potential first appear in ffcEF '■^^ order of 
|M^|C//A. Since it is enough for us to keep the condition of 
the weak CEF, |W^|?7/A should be smaller than the relevant 
energy scale of H^^^, i.e., the Hund's rule interaction in the 
j-j coupling scheme. Thus, we obtain another realistic condi- 
tion of \W\U/X<^E2. Here we note that the right-hand side 
is E2, not U. On the other hand, in the left-hand side, U is 
needed, since aU the Coulomb interaction terms contribute to 
eq. (3 1) in the intermediate processes. From eqs. (9) and (23), 
we obtain the revised condition of |Ty | / A<0.04 in the present 
parameter choice. In comparison with one of the original con- 
ditions of the weak CEF (|VF |/A<Cl), the range of A in which 
Hes works seems to be narrow, but when we compare the con- 
dition of |W|/A<0.04 with that of C//A<1, we understand 
that the range of A becomes wide so as to include the realistic 
parameter space. For instance, for the case of | VK|=10^''eV, it 
is allowed to take A in the order of 0. leV, which is the reahstic 
value for rare-earth ion. 

2.4 Stevens Hamiltonian in the weak CEF 

In the previous subsection, we have set the effective model 
-ffeff due to the expansion in terms of 1/A from the limit 
of A=oo. In order to assess the applicabihty of the effective 
model from the quantitative viewpoint, it is necessary to com- 
pare the results of HcS with those of H. For the purpose, in 
addition to the direct comparison between H and iJoS , it is 
useful to introduce the CEF Hamiltonian for Oh symmetry in 
the region of weak CEF, since we always consider the case 
of Z7>|iy| and A>|VK|. The CEF Hamiltonian is conven- 
tionally expressed by using the method of Stevens' operator 
equivalent as' 

Hs = Bl{n,J){dl + hdt)+Bl{n,J){6l - 210^), (32) 

where B^{n, J) and 0| denote, respectively, the CEF param- 
eter and the Stevens' operator equivalent for n and J. We call 
Hs the Stevens Hamiltonian. The matrix elements of for 
any value of J have been tabulated by Hutchings.^ Here we 
explicitly show the values of n and J in the parentheses of 
the CEF parameter, in order to distinguish them from B^ and 
B^ for J=£=3 in eq. (6). Note that Hs is the effective Hamil- 
tonian for the multiplet specified by J for any values of U 
and A, as long as they are much larger than \ W\. In fact, we 
have checked that the energy levels of H are correctly repro- 
duced by Hs with satisfactory precision for {/(::^|VK|) and 
A(>|M^|), except for the case of n=3 and A=oo. 

Effects of U and A appear in the CEF parameters, B^ {n, J) 
and Bg(n, J), which are expressed by 

Biin, J) = A4r^)(i^f\ B°{n, J) = Ae{r^Wf\ (33) 

where is the parameter depending on materials, and 
7j"^ are the so-called Stevens factors, which are coefficients 
appearing in the method of Stevens' operator equivalent,' 
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n 


J 


fc4(n, J) 


fc(i(n, J) 


L5 


j-j 


L5 


j-j 


1 


5/2 


11/7 


11/7 








2 


4 


-2/11 


-11/49 


-68/1155 





3 


9/2 


-340/4719 





1615/44044 





4 


4 


476/4719 


11/49 


-646/11011 





5 


5/2 


13/21 


-11/7 









Table I. Coefficients k4{n, J) and ke{n, J) for n=l~5 both in the LS 
and j-j coupUng schemes. 



and {r'') denotes the radial average of local /-electron wave- 
function. Note that in general, /? j"'' and 7j"'' depend on the 
Coulomb interaction and spin-orbit couphng, since they are 
determined by the nature of the ground-state multiplet spec- 
ified by J. In the present paper, the CEF potentials are al- 
ways given by and Bq. Thus, we express B^{n, J) and 
B^{n,J) as 

B°(n, J) = kiin, J)B2, Bg(n, J) = ke{n, J)Bl (34) 

By assuming that Ak and (r'^) are not changed in the same 
material group, we obtain 

k4{n, J) = /3S"VA. ke{n, J) = (35) 
where and 7^ for one / electron with 1=2) are given by^ 

= 2/(45 • 11), 7£ = -4/(9 • 13 • 33), (36) 

respectively. 

First let us consider the hmit of U=oo, i.e., the LS cou- 
pling scheme. We easily obtain fc4(n, J) and kQ{n, J) due to 
simple algebraic calculations by using the values of /Jj"^ and 
7j"'' for the LS coupling scheme.^ For n=l and J=5/2, from 
/^Vz =2/(45 -7), we easily obtain 

fc|'^(l,5/2) = ll/7, (37) 
as shown in eq. (21). For n=2 and J=4, we find 

fcf^(2,4)-/3f //3^--2/ll, .38. 
A:e^^(2,4)=7fV7£ = -68/1155. 

We can obtain k^^{n, J) and kg^{n, J) for n>3, as shown 
in Table 1. 

Next we consider another hmit of A=oo, i.e., the j-j cou- 
pling scheme. As long as we consider the situation of the 
weak CEF, it is possible to obtain the CEF energy levels by 
using the Stevens Hamiltonian Hs, even in the limit of X=oo, 
provided that B^{n, J) and B^{n^ J) are correctly evaluated. 
First we remark a couple of issues which can be understood 
without calculations: (i) Since B^ does not appear in the j-j 
coupling scheme due to the symmetry reason, we always ob- 
tain kg"-' {n, J)=0 for n<7 in the j-j coupling limit, (ii) For 
n=l, ki~^ (1, 5/2) is equal to fcf^(l, 5/2). 

Concerning ^4(77,, J) for n=2^5 in the j-j coupling limit, 
we evaluate the value in the following procedure: Let us 
consider the case of h=2 as a typical example. We note 
|2, 4, 4)=/^ g^2.fi 3/2!'^) '^he j-j coupling scheme, where 
|n, J, Jz) generally denotes the eigenstate determined by 
Hso+Hint, Jz is the z component of total angular momentum 
J, and |0) is the vacuum state. In the j-j coupling scheme, 
\n, J, Jz) for n<6 indicates the eigenstate of H^^l. Note that 
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the bra vector is defined as (n, J,Jz\- Then, we evaluate the 
CEF matrix element as 

(2, 4, 4|ijg|F|2, 4, 4) = -120i?°(l, 5/2) = -(1320/7)B°, 

(39) 

from eqs. (20) and (21). Since the same matrix element is 
given by 840i?4(2, 4) from the Hutchings table,^ we obtain 

fcr^(2,4) = -11/49, (40) 

in the j-j coupling limit for n=2 and J=4. Also for 
n=3~5, we can evaluate fc^^-* (n, J) in the similar way and 
the results are listed in Table 1. It should be noted that 

fci-^ (3,9/2)=0. Since |3, 9/2, 9/2)= fl,^jl,/jl,/,\0), 

we easily find (3, 9/2, 9 /2\H^^p\3, 9/2,9 /2)=0. Namely, for 
n=3 and J=9/2, the lowest-order contribution of the CEF po- 
tential is in the order of W'^, not in the order of W, in the j-j 
coupling scheme. This point will be discussed later again. 

For the intermediate coupling region in which both U and 
A are finite, we evaluate numerically k4{n, J) and ke{n, J) 
by deriving Hg from the original Hamiltonian H. First we 
diagonalize Hgo + -f^int as 

{H,o + Hint)\n, J, Jz) = E{n, J)|n, J, J,). (41) 

Note that the ground state has (2 J-|-l)-fold degeneracy, when 
the CEF potential is ignored. Then, we consider the CEF po- 
tential of Oh symmetry in the first order of W and the Hamil- 
tonian Hg is expressed in the matrix form as 

Hs{Jz, J'z) = {n, J, Jz\HcMn, J, J'z)- (42) 

Since the matrix elements of Hg have been already hsted for 

each value of ,/ by using B^in, J) and BQ{n, J),^ we can 
numerically obtain B^{n, J) and -Bg(n, J) from eq. (42) for 
any values of U and A for a given value of n. 

For comparison, we also evaluate B^{n, J) and BQ{n, J) 
from Hes. In this case, the procedure is almost the same as 
above. First we diagonalize H^s at W=0 as 

H^fi{W = 0)\n, J, Jz) = E^sin, J)\n, J, Jz). (43) 

Then, the Hamiltonian H^^ from ffgff for small W is ex- 
pressed in the matrix form as 

HfiJz, J'z) = {n, J, Jz\H,siW)\n, J, J'z) 

(44) 

- Ees[n, J)5j^j'^. 

In order to obtain B2{n, J) and Bl{n^J) numerically, we 
compare the matrix elements of eq. (44) with the table of 
Hutchings. 

3. Results 

3.1 p states 

First we explain the case of n=2 in detail. Before proceed- 
ing to the result of iJcfi for n=2, let us examine the results of 
the LS and j-j coupling schemes. Then, readers can under- 
stand how the j-j coupling scheme will be improved by i/off ■ 
In Figs. 2(a) and 2(b), we show the CEF energy levels ob- 
tained by the direct diagonalization of H for {U, A)=(10^,l) 
and (1,10"''), which correspond to the LS and the j-j coupling 
hmits, respectively. Note that the origin of the energy is ap- 
propriately shifted for convenience. 

In Fig. 2(a), the situation corresponds to the LS coupling 
scheme, in which first we obtain the ground-state level as 
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Fig. 2. (a) CEF energy levels in the coupling scheme forn=2. We diag- 
onalize H for f7=10^ and A=l. (b) CEF energy levels in the j-j coupling 
scheme for n=2. Here we set (7=1 and A=10^ for H. (c) k^{2, 4) and (d) 
fc6(2, 4) as functions of A/f7 for n=2. Broken curves denote A;4(2, 4) and 
fc6(2,4) evaluated from //eft- (e) Overlap integral as a func- 

tion of XjV for 71=2. Inset indicates the result in the intermediate coupling 
region in a magnified scale. The overlap is evaluated by the diagonaliza- 
tion of Hso + //int ■ (f) CEF energy levels of //eft (black solid curves) and 
those of H (red broken curves) for ?i=2, (7=1, and A=0.2. 



with S=\ and L=5 due to the Hund's rules in / shells. Upon 
further including the spin-orbit interaction, the ground state 
is specified by J=4 expressed as ^Ri^ in the traditional nota- 
tion. After considering the CEF potential, we find nine eigen 
states, singlet, V\ doublet, and two kinds of triplets (r4^ 
and Fg").^ It is found that the eigen energies of i7s for n=2 
and J=4 agree perfectly with those in Fig. 2(a), when we use 
yfcf ^(2, 4) and k^^{2, 4) in Table I. 

In Fig. 2(b), we show the results of H in the limit of the j-j 
coupling scheme. We have checked that these CEF energies 
are the same as those of H^a^ for n=2. Note that the results 
seem to be different from those in Ref. 4 at the first glance, but 
it is simply due to the negative sign of W . As mentioned in the 
previous section, even for the j-j coupling scheme, as long as 
we consider the situation of the weak CEF, it is possible to ob- 
tain the CEF energy levels by using the Stevens Hamiltonian 
iJs- We obtain the same results as those in Fig. 2(b) by diag- 
onalizing iJg with 5^(2, 4)=-(11/49)B2 and Bg(2,4)=0. 

The energy levels in Fig. 2(b) can be understood more in- 
tuitively as follows. In the limit of A=oo, the effective model 
iJcff is simply reduced to H^a^ , in which the only relevant 
CEF parameter is B\, leading to the splitting between Fy and 
Fg" levels for n=\. Thus, we obtain the eigen states for n=2 
by accommodating a couple of electrons in the potential for 
71=1 . For x<0 with positive B\, we find the V\ singlet ground 
state, including doubly occupied Ff orbitals, since Fy dou- 



blet is the ground state for n=\ in the region of x<0. The first 
excited state is F^ triplet, which is formed by Fy and Fg^ 
electrons. For a;>0, since the ground state for n=l is the Fg^ 
quartet, the ground state for n=2 is F^ triplet composed of a 
couple of F^ electrons, stabilized by the Hund's rule interac- 
tion. The first excited state in this region is F;^ doublet, which 
contains the component of a couple of Fg electrons. 

It is true that all the states appearing in the LS coupling 
scheme can be also found in the j-j coupling scheme. In par- 
ticular, at x=±l, the CEF energy levels agree well with those 
in Fig. 2(a). However, when we compare Figs. 2(a) and 2(b) 
in the region of I^It^I, the j-j coupling results seem to be dif- 
ferent from those of the LS coupling scheme, since the effect 
of Bg is not included, as already mentioned in the introduc- 
tion and in Ref. 4. This is quite natural, if we recall the fact 
that there exist no contributions of -Bg in the Hutchings table 
for J=5/2. This point is improved when we consider the 1/A 
corrections, as shown later 

In Figs. 2(c) and 2(d), we show fc4(2,4) and fcg(2,4), re- 
spectively, as functions of Xjl] . Since the results have been 
found to depend only on the ratio of A/[/, we change A and 
U appropriately by keeping the condition of the weak CEF, 
i.e., A>|VK| and C/>|M^|. In actual calculations, for A/C/>1, 
we gradually increase A from unity by setting [/=!, while 
for A/C/<1, U is increased from unity for the fixed value 
of A=l. Both in the limits of \/U=0 (LS coupling scheme) 
and \/U=oo (j-j coupling scheme), the numerical results 
correctly approach the analytic values shown in Table I. For 
the intermediate coupling region of X/U^Q.l, fc4(2,4) and 
fc6(2,4) smoothly change between the LS and j-j coupling 
values. We note that fc4(2,4)<0 and fc6(2,4)<0, suggesting 
that the sign of (2, 4) and Bg (2, 4) should be different from 
that of B^. In order to set 5^(2, 4)>0 and B°(2, 4)>0 for the 
comparison with the results in Ref. 2, the sign of W is taken 
as negative. 

Here we remark that fc4(2,4) and fc6(2,4) become good 
measures to show how the situation is close to the LS or the 
j-j coupling scheme. For actual materials, it is considered that 
U is about leV, while A is in the order of 0.1 eV, leading to 
X/U in the order of 0.1. From Figs. 2(c) and 2(d), ^4(2, 4) 
and fee (2, 4) in such a region are close to neither the value of 
the LS nor that of the j-j coupling scheme. Namely, the ac- 
tual situation corresponds to the intermediate coupling region. 
In particular, we note that the LS coupling scheme is differ- 
ent from the actual situation, in contrast to the naive expec- 
tation for the validity of the LS coupling scheme. One may 
consider that we can fit the results of the intermediate cou- 
pling region by changing B\{n, J) and B'^{n, J) in the LS 
coupling scheme. It is true that we can change A4 and Aq as 
fitting parameters, but it is not allowed to fit the Stevens fac- 
tors, which are determined by U and A under the symmetry 
constraint. Thus, the actual situation of the intermediate cou- 
pling region is different both from the LS and j-j coupling 
schemes. 

However, the above results on the smooth change of 
-B2(2,4) and i?g(2,4) seem to indicate that the LS and j-j 
coupling schemes are continuously connected. In fact, when 
we evaluate the overlap integral concerning the ground-state 
wavefunction, the magnitude is continuously changed be- 
tween the LS and j-j coupling schemes. In Fig. 2(e), we 
show the overlap integral as a function of X/U. 
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When we decrease \/U, we find that the overlap integral is 
gradually decreased, but even at the limit of U=oo, we still 
find {^Ls\^j-j)=0M5. Namely, the /^-state composed of 
a couple of electrons in the j=5/2 sextet becomes the good 
approximation for the LS coupling scheme, even though the 
j=7/2 octet is simply discarded. 

In principle, it is possible to approach the intermediate cou- 
pling region either from the LS or the j-j coupling limit, as 
mentioned in Sec. 2.2. Depending on the nature of the prob- 
lem, we can prefer one of them, but in this paper, the j-j 
coupling scheme is chosen. Here we remark that 
is close to unity in comparison with in the in- 

termediate coupling region. In fact, for X/U=0.l, we find 
(*|*^_^)=0.972, while (*|*ls)=0.957. 

Now we discuss the results of Hcs- In Figs. 2(c) and 2(d), 
broken curves denote ^4(2,4) and fee (2, 4) evaluated from 
eq. (44), respectively. We clearly observe that the results are 
actually improved from the values of the j-j coupling limit. 
In addition, it is found that they are close to the solid curves 
even for X/U in the order of 0. 1, i.e., in the intermediate cou- 
pling region. This fact supports the previous statement that 
-ffoff works even for A in the order of 0. 1 eV, when we take U 
in the order of eV. 

In Fig. 2(f), we show the energy levels of HcS by black 
solid curves for A=0.2 eV and U=l eV. For comparison, we 
show the results of H by red broken curves for the same pa- 
rameters. Note that we directly diagonalize Hcs and H to de- 
pict solid and broken curves. In the first impression, the re- 
sults of HcS are similar to those of H. Since the effect of Bq 
is now included efficiently in the effective interaction eq. (31), 
the CEF energy levels of H are well reproduced. In particu- 
lar, we can obtain Fj^ doublet ground state, which does not 
appear in the limit of A=cx3. 

Although we do not have perfect agreements between the 
results of iJcff and H in Fig. 2(f), it should be noted that we 
do not adjust any parameters. In this sense, it is concluded that 
the characteristic features of the CEF energy levels of H are 
well captured by Hcs- Readers may consider that it is possible 
to fit the results of H by the LS coupling scheme due to the 
adjustment of the CEF parameters. However, the adjustable 
parameters, A4 and in eqs. (33), are considered to depend 
on materials, not on U and A. As remarked above, we cannot 
adjust the ratios of the Stevens factors, /c4(n, J) and kQ{n, J), 
since they are determined by the symmetry requirement and 
the values of U and A. Namely, the difference in the inter- 
mediate coupling region between the original Hamiltonian H 
and the LS coupling scheme cannot be improved essentially, 
as long as we consider the limit of U=oo. If we actually in- 
tend to improve the LS coupling scheme, it is necessary to 
consider, for instance, the expansion in terms of 1/U from 
the limit of U=oo. This is an alternative way to construct the 
effective model in the intermediate coupling region. 

3.2 states 

Next we consider the case of n=3 by following the same 
discussion flow as in the previous subsection for n=2. In 
Figs. 3(a) and 3(b), we show the CEF energy levels obtained 
by the direct diagonalization of H for (U, A)=(10^, 1) and 
(1,10^), respectively. In the LS coupling scheme, the ground 
state multiplet is characterized by ^19/2 iL=6, 5=3/2, J=9/2). 
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Fig. 3. CEF energy levels (a) in the LS coupling and (b) the j-j coupling 
schemes for n=3. To express virtually the LS and j-j coupling schemes, 
we set (U, A)=(105, 1) and (1, 10^) for H in (a) and (b), respectively, 
(c) fc4(3, 9/2) and (d) ^6(3, 9/2) as functions of \/U for n=3. Broken 
curves denote ^4(8, 9/2) and fcgCS, 9/2) evaluated from H^g. (e) Over- 
lap integral ('I'I'I'j-j) as functions of \/U for 7i=3. Inset indicates the 
result in the intermediate coupling region in a magnified scale, (f) CEF 
energy levels of H^g (black solid curves) and H (red broken curves) for 
n=3, U=\, and A=0.2. 



The 10-fold degenerate state is split into two and one F^ 
states under the Cubic CEF.^ On the other hand, in the limit of 
the j-j coupling scheme, the situation seems quite different. 
As shown in Fig. 3(b), the Fg quartet always becomes the 
ground state, while the energy for the 6-fold degenerate ex- 
cited states including another quartet and Fg^ doublet do 
not depend on the CEF potential. We find that the curve for the 
ground state quartet is quadratic in terms of the CEF potential, 
suggesting that the contribution in the order of W vanishes 
in the limit of the j-j coupling scheme, as suggested from 
fcr^ (3, 9/2)=fc^~^ (3, 9/2)=0. Namely, the Stevens Hamilto- 
nian i7g does not work in this case. However, it is found that 
the results of H^i^ for n=3 correctly reproduce the CEF en- 
ergy levels in Fig. 3(b). This fact clearly indicates that it is 
necessary to calculate the energy levels by diagonalizing si- 
multaneously the Coulomb interaction and the CEF potential 
terms in the j-j coupling scheme. As already mentioned, we 
always perform the direct diagonalization of HcS in the j-j 
coupling scheme. The Stevens Hamiltonian in the j-j cou- 
pling scheme is used only for the purpose to estimate {n, J) 
and kQ{n, J). 

Figures 3(c) and 3(d) denote ^4(8, 9/2) and ke{3, 9/2), re- 
spectively. We find that the values at the limits of U=oo and 
A=oo are exactly the same as the analytic values in the LS and 
j-j coupling schemes (see Table. I). In particular, we con- 
firm that fc^~^(3,9/2)=A:^~-''(3,9/2)=0, consistent with the 
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Fig. 4. (a) CEF energy levels in the LS coupling scheme for n=4, obtained 
by diagonalizing H for (7=10^ and A=l. (b) CEF energy levels in the 
j-j coupling scheme for n=A. Here we set U=\ and A=10^ for H. (c) 
fc4(4,4) and (d) fcg(4,4) as functions of \/U for n=4. Broken curves 
denote A;4(4, 4) and fc6(4, 4) evaluated from i^ef^■ (2) Overlap integral 
(^I'J'j—j) as functions of \/U for n=4. Inset indicates the result in the 
intermediate coupHng region in a magnified scale. Inset shows the result 
in the intermediate coupling region, (f) CEF energy levels of H^g (black 
solid curves) and H (red broken curves) for n=4, U=\, and A=0.2. 



quadratic behavior of the CEF energy in Fig. 3(b). In Fig. 3(e), 
we depict the overlap integral as a function of \/U . When 
we decrease \/U, the overlap integral is decreased and at 
the limit of U=oo, we find (\l'L5|*j-j)=0.674. It seems to 
be larger than we have naively expected. In the intermediate 
coupling region, we find (^'|^'j_j)=0.862 for A/[/=0.1 and 
(*|4'j_j)=0.932 for A/C/=0.2. These results indicate that it 
is meaningful to construct the /■^-electron state only by using 
the j=5/2 sextet in the intermediate coupling region for n=3. 

In Figs. 3(c) and 3(d), we also show fc4(3,9/2) and 
fc6(3,9/2) evaluated from eq. (44) by broken curves. They 
deviate from the solid curves around at A^l, but as shown in 
Fig. 3(f), even for A=0.2 and U=\, characteristic features of 
the CEF energy levels of H are reproduced by H^s- Thus, in 
the same approximation level when we use the LS coupling 
scheme, we can exploit H^s for the purpose to fit the exper- 
imental results by adjusting and in eqs. (33). There is 
an advantage of Hcs that the common CEF parameters can be 
used, even when the local /-electron number is changed. It is 
concluded that our modified j-j coupling scheme works well 
to consider the multi-/-electron state. 

3.3 states 

Encouraged by the results of n=2 and 3, let us further pro- 
ceed to the case of n=4. In Figs. 4(a) and 4(b), we show the 
CEF energy levels for n=A obtained by the direct diagonal- 
ization of H for {U, A)=(10^, 1) and (1,10^), respectively. For 



n=4, the ground state multiplet is characterized by ^1^ (L=6, 
S'=2, J=A). The nontet is split into singlet, doublet, 
and two kinds of triplets (FJ and V^)} Since the value of J 
for n=4 is equal to that of n=2, the CEF states are in com- 
mon with Fig. 2(a), but the x dependence is different, since 
the sign of /cf '^(4, 4) is different from fc|"^(2, 4). Due to the 
same reason, we find that the x dependence of the CEF energy 
levels in Fig. 4(b) is reversed from that of Fig. 2(b) in the j-j 
coupling scheme. 

In Figs. 4(c) and 4(d), we show k^^AA) and fc6(4,4) by 
solid curves, indicating that the values in the LS and the j- 
j coupling limits are correctly reproduced with the smooth 
changes between them. When we compare them with the re- 
sults for fc4(4, 4) and fc6(4, 4) evaluated from eq. (44) (broken 
curves), they begin to deviate from the solid curves around 
at A=2^3. As shown in Fig. 4(e), at the limit of [/=oo, 
we obtain (^i5|vl'j_j)=0.465 for n=4. In the intermediate 
coupling region, we find (5'|5'j_j)=0.692 for X/U=Q.l and 
(*|*j_j)=0.836 for X/U=0.2. Even for n=4, in the interme- 
diate coupling scheme, the ground state is well approximated 
by the j-j coupling scheme. Then, it is still possible to repro- 
duce the results in the intermediate coupling region. In fact, 
as shown in Fig. 4(f), the CEF energy levels of H look similar 
to those of Hes for U=l and A=0.2. 

3.4 p states 

Now we move on to the case of n=5. In Figs. 5(a) and 5(b), 
we depict the CEF energy levels due to the direct diagonaliza- 
tion of H with n=5 for (U, A)=(10^ 1) and (1,10^^), respec- 
tively. For n=5, the ground state multiplet is characterized by 
^H5/2 (L=3, 5=5/2, J=5/2). The sextet is split into Ff dou- 
blet and Fg^ quartet.^ 

Here readers may consider that the CEF energy levels for 
n=5 are simply the same as those in the case of n=l, but we 
should not simply conclude it. Namely, as easily understood 
from Figs. 5(a) and 5(b), the x dependence of the CEF en- 
ergy levels is just reversed between the LS and j-j coupling 
schemes. It should be noted that this phenomenon is not due 
to the approximation, but the intrinsic feature of the origi- 
nal Hamiltonian H. Note also that for n=l and 5, the only 
relevant CEF parameter is i?4(n, J), since the sixth-order 
CEF potential terms do not appear in the space of J =5 12. 
In the limit of the LS coupling scheme, since fc4(5,5/2) 
and fc4(l, 5/2) have the same sign, the ground state charac- 
ter should not be changed between the cases of n=\ and 5. 

On the other hand, in the limit of the j-j coupling scheme, 
the X dependence of the CEF ground state of n=5 is reversed 
from that of n=l, since the sign of fc4(5, 5/2) is different from 
that of fc4(l, 5/2) in the j-j coupling limit, as shown in Ta- 
ble I. We can understand the reason more intuitively from the 
electron-hole relation on the basis of the j-j coupling scheme. 
Namely, to obtain the /''-electron state in the j-j coupling 
scheme, we accommodate five electrons in the one-electron 
potential levels, for instance, such as Fy ground and Fg ex- 
cited states for ^^(l, 5/2)>0 (x<0 and W<0). Then, we ob- 
tain ground and F^ excited states for the case of ?t,=5 
in the same region of x. In Fig. 5(c), we depict fc4(5,5/2) 
as a function of X/U for n=5. We confirm that fc4(5,5/2) 
smoothly changes from —11/7 at A=oo to 13/21 at U=oo. 

When we compare the solid curves with the broken ones 
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Fig. 5. (a) CEF energy levels in the LS coupling scheme for n=5. We di- 
agonalize H for U=10^ and A=l. (b) CEF energy levels of H in the j-j 
coupling scheme for n=5. Here we set U=l and A=10^. (c) fc4(5,5/2) 
as a function of X/U for n=5. Broken curve denotes the result estimated 
from H^f[. (d) Overlap integral (^I'l/j-j) as a function of X/U for n=5. 
Inset denotes the result in the intermediate coupling region in a magnified 
scale. CEF energy levels of H^f{ and H for (e) U=l and (f) [7=5.5 with 
n=5 and A=0.2. In both figures, we depict the results of i?cft and H by 
black solid and red broken curves, respectively. 



for ^4(5, 5/2) estimated from eq. (44), the deviation begins 
to occur at A between 1 and 10. Namely, the degree of ap- 
proximation seems to be worse in comparison with the case 
of n=4. In fact, we obtain the smaller overlap integral for 
n=5, (^'rsl^' 



=0.280, as shown in Fig. 5(d). However, 
in the intermediate coupling region, we find (5'|5'j_j)=0.622 
at A=0.1 and (*|«'j_/)=0.857 at A=0.2, for a fixed value of 
U=l. At A=0.2, {^l^j^j) for n=5 is larger than that for n=4, 
as discussed in Sec. 2.2. In any case, the ground-state wave- 
function in the intermediate coupling region is still approx- 
imated well by the states constructed from the j-j coupling 
scheme for n=5. We can reproduce the results of H by -ffeff 
for U=l and U=5 at A=0.2, as shown in Figs. 5(e) and 5(f). 
In particular, in the effective model of Hctf, we can reproduce 
the interchange of the ground states, when we change the ratio 
of X/U. Note, however, that we cannot quantitatively repro- 
duce the critical value of U at which the ground states are 
interchanged. 

In order to understand the interchange of the ground state 
for the case of n=5 within the LS coupling scheme, the only 
way is to change the sign of W phenomenologically, since 
k^^{5, 5/2) is fixed due to the symmetry requirement. How- 
ever, in our modified j-j coupling scheme, we can consider 
the microscopic origin of the change of the CEF parame- 
ter due to the competition between the Coulomb interaction 
and the spin-orbit coupling. It is an advantage of our effec- 
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Fig. 6. CEF energy levels of H^g (black solid curves) for (a) n=l and (b) 
n=2 for Tjj symmetry with j/=0.3. The spin-orbit interaction is set as (a) 
A=0. 1 and (b) A=0.2. In the panel (b), we also show the results of H by 
red broken curves for the same parameters. 



live model, in addition to the reproduction of the CEF energy 
levels of H. 

As for the change of the sign in the CEF parameter, when 
we consider the different material groups, we may choose 
positive or negative W depending on the kind of materials. 
However, for the same material group with the same rare- 
earth ion, it is difficult to imagine that the sign of the CEF 
parameter is changed, although the magnitude may be dif- 
ferent due to the change of ligand ions. Rather, as found 
in the present calculations, it seems natural to understand 
that the CEF ground state is converted due to the effect of 
the Coulomb interaction and/or the spin-orbit coupling. This 
is considered to explain a possible conversion of the CEF 
ground state in Sm-based filled skutterudites, as mentioned 
in the next subsection. 

3.5 Th symmetry 

In order to show the effectiveness of our modified j-j cou- 
pling scheme, let us also discuss the CEF energy levels for 
filled skutterudites with Th symmetry.^'^'^^ For the purpose, it 

i2 
6 



is necessary to add extra terms^^ in Bj^.m' as^ 



^3,1 = B_3,-i = 24VT5BI, 



B2.0 - 
Bi,-i 



B-2.0 - 
= -B3. 



48^15^6^ 
= 360B|. 



(45) 



By following Ref. 22, we express Bq as i?|=VFy/F*(6) with 
-F'*(6)=24 for J=3. Here we set y=0.3 as a typical value for 
filled skutterudites. 

In Figs. 6(a) and 6(b), we show the energy levels for n=l 
and 2 by diagonahzing H^tf for A=0.1 and 0.2, respectively, 
with [7=1 . For the case of n=2, we also depict the results of H 
by red broken curves. For n=l, we have doublet and Fg^ 
quartet states, which are essentially the same as doublet 
and quartet for Oh symmetry. For n=2, we find remark- 
able difference from the case of Oh symmetry. As already 
mentioned in Ref. 22, two triplets and in Oh symme- 



try are mixed and they are r^*-^-* and T^^^"' in Th symmetry. 
It is found that the results of H in the intermediate coupling 
region are well reproduced by H^s. It is quite natural, since 
the effect of Bg is effectively included in H^s. It is empha- 
sized here that such a characteristic issue of Th symmetry is 
correctly reproduced in our modified j-j coupling scheme. 

For n=5, we have (T^ in Oh) doublet and Fgy (F^ 
in Oh) quartet states. Thus, the results are essentially the 
same as those in Figs. 5(e) and 5(f). As mentioned above, 
for Sm-based filled skutterudites, it has been recently pointed 
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Fig. 7. CEF energy levels of J?^"' (black solid curves) for (a) n=13 and 
(b) 71=12 with U=l and W=—10~'^. For comparison, in the panel (b), we 
also show the results of H in the LS coupling scheme (C/=10^ and A=l) 
by red broken curves. 
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out a possibility that the ground states are changed between 
doublet and Fg^ quartet. For SmRu4Pi2, SmOs4Pi2, and 
SmOs4Sbi2, it has been considered that Fg^ quartet is the 
ground state. ^'^"^^ On the other hand, for SniFe4Pi2, a pos- 
sibility of FjT doublet ground state has been suggested.^''-'"' 
It can be interpreted that the conversion of the CEF ground 
state is due to the change of the Coulomb interaction and/or 
the spin-orbit coupling, which naturally leads to the change 
in the sign of B^. This point will be discussed in detail else- 
where with numerical results on the possible multipole state 
of Sm-based filled skutterudites.^'* 

4. Discussion and Summary 

In this paper, we have proposed the effective model to con- 
sider systematically the multi-/-electron states by using the 
expansion in terms of 1/A. We have shown that the results 
in the intermediate coupling region for n=2~5 can be repro- 
duced well by the effective model and the appUcability of the 
model has been clarified. The CEF state for Th symmetry can 
be also reproduced. 

We have stated that Hi^"^ cannot reproduce even qualita- 
tively the energy levels of the intermediate coupling region as 
well as the LS coupling scheme, since the effect of Bg is not 
included. Here readers may have a question: If the effect of 
Bq is included in the one-electron potential, does the zeroth- 
order term mimic the LS coupling results? The answer is yes. 
This can be clarified by considering the cases of n=13 and 12, 
i.e., the situations of one and two holes in the j=l/2 octet. For 
j=ll2, the zeroth-order term of 1/ A is given by 



(0) _ f,b,b A f 



Ejbb,bb ft ft f f 

Pi P2 -,0304 J bPi J b/32 ■' ^/53 ■' ' 

/3i~/34 



(46) 



where iJ^ « is given by the CEF potential for j=7/2 and 



rbb,bb 



'Pi, 132 



^/5i/?2 0304 Coulomb interaction among j=7/2 octet. 

In Figs. 7(a) and 7(b), we show the energy levels of by 
black solid curves for ni,=l and 6, corresponding to n=13 and 
12, respectively, where ni,=n — 6 and ni, denotes the electron 
number in the j=7/2 octet. Note that the results are obtained 
by the diagonalization of -ff^"^ . For comparison, for the case 
of n=l2, we also show the results in the LS coupling scheme, 
by setting C/=10^ and A=l for H. In the j=7/2 octet, the ef- 
fect of Bg is already included in the one-hole case (n=13). 



as easily understood from the Hutchings table for J =112. For 
the two-hole case (n=12), as shown in Fig. 7(b), even with- 
out considering the 1/A-correction term the CEF en- 
ergy levels of H^^ (black solid curves) are quite similar to 
those in the LS coupling scheme for J=6. Of course, in or- 
der to reproduce quantitatively the CEF energy levels in the 
realistic intermediate coupling region, it is necessary to in- 
clude the terms in the order of 1/ A for the j=7/2 octet, as has 
done in this paper for the j=5/2 sextet. However, for practical 
purposes such as the fitting of experimental results, it seems 
enough to use H^^'' for Tm and Yb compounds. When we 
further proceed to the cases of n<12, it is recommended to 
improve i/^^-* by considering the 1/A corrections. It is one 
of future issues, when we attempt to develop a microscopic 
theory of magnetism in heavy rare-earth compounds. 

As for a possible application of our modified j-j coupling 
scheme, we consider a direction to improve effectively the 
band-structure calculations for /-electron materials. In the 
relativistic band-structure calculations, all the one-electron 
CEF potentials for j=5/2 and 7/2 states are correctly in- 
cluded,^^ but the multi-/-electron state is not completely re- 
produced in comparison with the CEF levels of the LS cou- 
pling scheme. In order to construct the multi-/-electron state 
due to the one-electron basis, it is necessary to treat correctly 
the competition among Coulomb interactions, spin-orbit cou- 
pling, and CEF potentials, but in the band- structure calcula- 
tions, the effect of Coulomb interactions, in particular, the 
Hund's-rule coupling, is considered only partly within the 
mean-field approximation. An effective way to improve such 
a situation is to include the two-body CEF potentials dis- 
cussed here in the band-structure calculations with due care. 
It is one of future problems to develop a systematic way for 
the inclusion. 

Another issue is the CEF effect on exotic itinerant mag- 
netism and unconventional superconductivity of /-electron 
materials from the band picture. We can construct the ap- 
propriate many-body Hamiltonian for /-electron systems, by 
adding the hybridization term between / and conduction elec- 
trons or the hopping term of / electrons to the local /-electron 
term Hcff in the modified j-j coupling scheme. The parame- 
ters in the kinetic term are determined so as to reproduce the 
energy bands near the Fermi level. It is another future problem 
to analyze such a model by using numerical and/or analytical 
techniques. 

In summary, we have discussed the /"-electron states with 
n>2 on the basis of the effective model obtained by including 
the corrections in the order of 1/ A in the j-j coupling scheme. 
For n=2r-^5, the results in the realistic intermediate coupling 
region have been quantitatively reproduced in our modified j- 
j coupling scheme. By using HcS, we have also reproduced 
correctly the CEF energy levels for Th symmetry. In conclu- 
sion, for the consideration of multi- /-electron states, we can 
use the j-j coupling scheme with appropriate corrections in 
terms of 1/A. 
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